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Abstract 

In this paper we construct the trajectory-coherent states of a damped har- 
monic oscihator. We investigate the properties of this states. 
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I. INTRODUCTION 



At present the coherent states are widely used to describe many fields of theoretical 
physics [|l|. The interest and activity in coherent states were revived by the paper of Glauber 
0, who showed that the coherent states could be successfully used for problem of quantum 
optics. Recently Nieto and Simmons have constructed coherent states for particles in 
general potential, Hartley and Ray Q have obtained coherent states for time-dependent 
harmonic oscillator on the basis of Lewis-Risenfeld theory; Yeon,Um and George |Q have 
constructed exact coherent states for the damped harmonic oscillator. 

Some time ago Bagrov,Belov and Ternov |^ have constructed approximate (for h 
) solutions of the Schrodinger equation for particles in general potentials, such that the 
coordinate and momentum quantum - mechanical averages were exact solutions of the cor- 
responding classical Hamiltonian equations; these states were called trajectory- coherent 
(TCS). The basis of this construction is the complex WKB method by V.P.Maslov 

The aim of this work is to construct the trajectory-coherent (TCS) states of a damped 
harmonic oscillator by using the Caldirola - Kanai Hamiltonian and the complex WKB 
method. It is shown that this states satisfy the Schrodinger equation exactly. 

II. THE CONSTRUCTION OF TCS. 

Consider the Schrodinger equation 

indt"^ = (1) 

where the symbol of operator H - the function H{x,p,t) is arbitrary real and analytical 
function of coordinate and momentum. The method of construction TCS in this case has 
been described in detail in Q, hence we shall illustrate some moments only. For constructing 
the TCS of the Schrodinger equation it is necessary to solve the classical Hamiltonian system 

x{t) = dpH{x,p,t), p{t) = -d^H{x,p,t), (2) 

and the system in variations (this is the linearization of the Hamiltonian system in the 
neighbourhood of the trajectory x{t),p{t)) 

wit) = -H,p{t)w{t) - H,.{t)z{t), w{0) = 6, (3) 
zit) = H„{t)w{t) + HpS)z{t), z{{}) = 1, 

where H{x,p,t) is the classical Hamiltonian; 

Hxpif) = dxdpH{x,p,t) 

Hpx{t) = dpdxH(x,p,t) \x=x{t),p=p(t), 

Hxx{t) = dl^H{x,p,t) \x=x{t),p=p{t), 

Hppit) = dppH{x,p,t) \x=x(t),p=p(t), 

b is complex number obeying the condition Im6 > 0, and x{t), p{t) are the solutions of 
system (2). 
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Consider the damped harmonic oscillator []5|, [l0| , [Tl 



H{x,p,t) = exp{—'~ft){2m) ^p"^ + -exp{'~ft)mujQx'^ . (4) 



The Lagrangian and mechanical energy are given by 

.1 .o 1 



L = exp(7t)(-mx - -mc^gX ), (5) 
E = exp{—2'~ft){2m)^^p'^ + -mc^QX^. 



We first define the Cauchy problem for equation (1) 

I >|t=o= "^oix, t,h) \t=o= N exp{ih'^{po{x - xo) (6) 

+ -(x-xo)^)}, 

where Xq = x{t) \t=o, po = p{t) \t=Q . 

The function of WKB - solution type (TCS) § 

I >= ^o(a;, t, h) = iV$(t) exp{i^-^^(x, t)}, (7) 

where = {Imb{nh)-^y/\ = {z{t))-^/^, 

t 

S{x,t) = J {x{t)p{t) - H{x{t),p{t),t)}dt 



+p{t){x - x{t)) + ]^w{t)z-^{t){x - x{t) f, 

and the phase S{x, t) is the complex - valued function (ImS' > 0) is the approximate solution 
of the Cauchy problem (6) for the Schrodinger equation (1). We should note that in the 
case of the quadratic systems, for example, for Hamiltonian (4) the function (7) is the exact 
solution of the equation (1). 

Solving the differential equations (2), (3) we obtain 

x(t) = exp( — 7t) (2posinLi;t + m(2ci;coscc;t + 7sinu;t)a;o) , (8) 

2muj 2 

11/ 

p(t) = exp(-7t) ((27sincijt — 4u;cosu;t)po + ^(7^ + 4u;^)a:o sincut 

4a; 2 ^ 

w{t) = —- — exp(-7t) ((267 + m(7^ + 41^;^)) sincut — Abu cosut) , 

i f (2b + mj) smu!t\ 
z{t) = eM-^lt) cosu;t+ ^— , 
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It is easy to check that the function (7), where x{t),p{t), w{t), z{t) are defined by (8) satisfy 
the equation (1) exactly as for u;^ > as < {u ^ iuo) . 

Further, we define "annihilation" operator d(t) and "creation" operator a^{t) as 0: 

a{t) = {2hIinb)-^/^{z{t){p - p{t)) - w{t){x - (9) 

a+(t) = {2hluihY^/'^{z*{t){p - pit)) - w*{t){x - x{t))}. 

It easy to check that the creation and annihilation operators satisfy the usual Bose permu- 
tation rule 

[a, a+] = 1, [a, a] = [a+, a+] = 0. (10) 
The complete orthonormal set of the trajectory-coherent states (TCS) are defined as |0: 

I n >= (n!)-i/2(a+)" | > . (11) 
It is not difficult to check the relations 

< n \ m >= Sn,m, 0. I >= 0, (12) 

a+ I n >= {n + 1)^/^ | + 1 >, a\n >= (n)^/^ \ n - 1 > . 

Using the usual methods ^ we obtain the expression for coherent states (CS) 

I a >= A{a) I >, (13) 

where A{a) is the unitary operator 

A{a) = exp(Q;a^ — a*d), (14) 

a is the complex number; d~^,a are defined by (9). 

Besides, at is follows from (12)-(14), the functions | a > are eigenstates of the opera- 
tor a with eigenvalue a, i.e. 

a I a >= a \ a > . (15) 



III. QUANTUM - MECHANICAL AVERAGES AND UNCERTAINTY 

RELATIONS. 

Further we shall find the expressions for quantum - mechanical averages 

< X >TCSi < ^ >TCS) < P >TCS) < P >TCS) < ^ ^CSy 



where, for example, 

< X >Tcs—< n \ X \ n >, < X >cs—< a \ X \ a > . 

Purther,we present the relations expressing the coordinate and momentum operators in 
terms of the operators a 

X = x{t) - i{2Inib{h)-Y^/^z{t)a+ - z*{t)a}, (16) 

p = p(t) - i{2lrah{h)-^)-^''^{w{t)a+ - w*{t)a}. 
Using (5), (9)- (16), we obtain for quantum - mechanical averages: 

< X >TCS^ X{t), < p >TCS= P{t), (17) 

< X >cs^ x{t) - i{2Inib{h)-Y^/^{oi*4t) - 0iz*{t)}, 

< p >^^= p(t) - i{2Iinb{h)-Y'^/^{oi*w{t) - aw*{t)}, 

< x" >^cs= x^t) + ^(n + ^)\z{t)\^, 

< f >Tcs^ p\t) + -^(n + \)\w{t)\\ 

imo 2 

< X^ >cs=< X >ls 

< f >cs^< P >ls 

< E >Tcs= cxp(-27t)(2m)"y(t) 

x{exp(-27t)(2m)-V(^)r + -^^okWD^ 

2 



< E >cs^ exp(-27t)(2m) ^ < p 

1 2 2 ^ 

+-mu;o <x>cs + 



2Imb 

x{exp(-27t)(2m)"V(^)r + ^^^okWI^}- 
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Now, by calculating the uncertainty in x and p in the TCS and CS one finds: 

{Ax)l.cs =< < ^ >TCsf >TCS (18) 

{^P)lcS =< iP- < P >TCs)^ >TCS 

= h{n+ ^' ^ ^' 



2' Imb 



{Ax)l, =< {x- < X >cs)' >as= ^^1^, (19) 



2 Im6 



{Ap)l, =< {p- < p >cs)^ > 

So, the Heisenberg's uncertainty relations is expressed as 



h \w{t)\' 



1 2KM)I' 

2^ {Imby 



{Ax)l,,{Ap)l,, = h\n + ^)^ '7;iX^' , (20) 



(A^)L(AP)L - f ^ff^. (21) 

For the damped harmonic oscillator (4) we choose Re6 = (it is necessary for mini- 
mization of the uncertainty relations in the initial instant of time [0]), Imb = fimuj , where 
/i > shows the initial uncertainty of coordinate; and we denote 6 = '-f/2uj. 

By using (8),(16)-(21), in the case u"^ = Uq - > for (Ax)^^^, 
(Ap)Ls' (A£)cs' (Ap)cs we obtain 

(Ax)ypg = h{n -\ — ) exp{—'-ft){fimuj)^^{l + sin^u;t(^^ 

- 1) + ^sin2a;t}, 



(Ap)^^g = h{n H — ) exp(7t)/imc<j{l H — ^ sin^ ujt{26'^ 

+ 1 + n^e"^ - /i^) - ^ sin 2cut}, 



(Ai:)^g = -exp(-7t)(/imcj)"^{l + sin2cjt(6'^ 

- 1) + esi\i2ujt}, 

h 1 
(Ap)^^ = — exp(7t)/imt<;{l H — ^ sin^ co't(2^^^ 

+0^ + 1 + /i^e^ - /i^) - ^ sin 2cjt}, 
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and thus the uncertainty relations becomes 



{Ax)l,,{Ap)l,, = h\n + + g{t)l (22) 

{M)l,{Ap)l, = -{l + g{t)), 

where 

Q 

g{t) = {-{9^ + + I) sin^ ut (23) 
/i 

+— (^^ -/i^ + l)sin2wt}l 
2/i 

We should note,that in the special case ^ = 1 the formula (23) coincides with formula 
(14) from The function (23) is equal to zero,and,therefore, the minimization of the 
uncertainty relations has place in the instants of time 

ttA; 1 /i^ - 02 _ ^ 

tik = — , t2k = — arctan ■ 



u e{^j? + 9^ + 1) 

+ — ; k = 0,l,2,.... 

UJ 

The preceding equation can be solved for the parameter /x > in the case 

\9iBxiujt\ < 1. (24) 

and,therefore,by choosing parameter /i we can obtain the minimization for any instant of 
time t obeying the condition (24). 

In the case u"^ = — cJq > we obtain 

(^*)tcs = + ^) exp(-7t)(/im^)~^ 

X {1 + sinh^ ujt{9'^ + 1/ + 1) + 9 sinh 2ujt}, 

(^p)tcs = + ^) exp(7t)/imo; 

X {1 + ^ sinh' tutil - 29^ + 9^ + ^^9^ + ^i^) 
—9 sinh 2ujt], 

(Ax)^^ = —exY)[—^t){jj,mujy^{l + sinh'c<jt(0' 



+ 1) + 9 ^inh.2ujt} , 
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^exp(7t)Aimcj{l + 

+ ^^e"^ + l^'^)-0 sinh 2ut}, 



{Ap)ls = -exp{jt)fimuj{l + — sinh^a;i(l - 261^ 



and for the uncertainty relations we have 

iAx)l,,{Ap)l,,^h\n + lril + g{t)), 

{Ax)UAp)l,^^{l + g{t)), 

where 

e 



g(t) = {-{e^ + A*^ - 1) sinh^ out (25) 
A* 



Now the function g{t) (25) is equal to zero only for 

1 n'-d^ + l 

toi = , to2 = — arctanh- 



oj e{ij? + ^2 _ 1) • 

This equation can be solved for parameter fx in the case 

|^tanht,;t| < 1 (^ > 1), |^ tanhc;t| > 1 (0 < ^ < 1) (26) 

and we can obtain the minimization for any instant of time t obeying the condition (26) .We 
should notice, that the formula (23) can be easily obtained from (25) by means of replacement 
uj ^ iuj and, therefore ^ —'i-fJ', d —id- In a case of harmonic oscillator (7 = 0) as (23) 
and (25) bring to g{t) = 0, that coincides with results, received in Ref.l2. 
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